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Abstract 

We investigate quasi-Monte Carlo integration using higher order digi¬ 
tal nets in weighted Sobolev spaces of arbitrary fixed smoothness a € N, 
a > 2, defined over the s-dimensional unit cube. We prove that ran¬ 
domly digitally shifted order /3 digital nets can achieve the convergence of 
the root mean square worst-case error of order A““(log when 

/? > 2a. The exponent of the logarithmic term, i.e., (s — l)/2, is improved 
compared to the known result by Baldeaux and Dick, in which the expo¬ 
nent is sa/2. Our result implies the existence of a digitally shifted order 
/3 digital net achieving the convergence of the worst-case error of order 
A““(log which matches a lower bound on the convergence rate 

of the worst-case error for any cubature rule using N function evaluations 
and thus is best possible. 

Keywords: Quasi-Monte Carlo, numerical integration, higher order digital nets, 
Sobolev space 
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1 Introduction and the main result 

In this paper we investigate quasi-Monte Carlo (QMC) integration of functions 
defined over the s-dimensional unit cube. For an integrable function /: [0,1)'* —>■ 
R, we denote the true integral of / by 

/(/):=/' fix)dx. 
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QMC integration of / over a finite point set P C [0,1)® approximates /(/) by 

' ' X&P 

Here P is a mnltiset and so if an element occurs mnltiple times it is counted 
according to its multiplicity. The key ingredient for success of QMC integration 
is to construct good point sets depending on a function class to which / belongs. 
In the classical QMC theory, for instance, a class of functions with bounded 
variation in the sense of Hardy and Krause has been often considered [iniiii]. 
For this function class, the Koksma-Hlawka inequality states that the integration 
error is bounded by 


|/(/;P)-/(/)| <I/hk(/)P*(P), 

where Vhk(/) denotes the variation of / in the sense of Hardy and Krause, 
and D*{P) the star-discrepancy of P, see [TH Chapter 3]. This inequality 
motivates construction of low-discrepancy point sets. We refer to [71 Chapter 8 ] 
for several explicit constructions of point sets whose star-discrepancy is of order 
iV“^(logiV)®“^, where N denotes the number of points, i.e., N = |P|. 

One of the recent interest in the research community is to consider a function 
class consisting of smooth functions and to construct good point sets which 
achieve higher order convergence of the integration error in that function class, 
see for instance mniiiiiiin]. a function space of our particular interest in this 
paper is a weighted Sobolev space Pa ,7 of arbitrary fixed smoothness a for a set 
of non-negative real numbers 7 = (7u)uc{i,...,s}- Here a > 2 is a positive integer. 
(We shall give the precise definition of Pq,^ in Subsection 12.Il l This space has 
been studied, for instance, in [DI3I7] in the context of QMC integration. The 
breakthrough in this research direction was made by Dick and his collaborators 
[Dia E121312], who provide us with an explicit construction of good point 
sets called higher order digital nets achieving almost optimal convergence of the 
integration error of order A^““(log for some c(s,a) > 0. (We shall give 

the definition of higher order digital nets in Subsection 12.21 1 The above order 
of convergence a is best possible up to some power of a log N factor. 

A thorough analysis on the exponent c(s, a) has been recently done for peri¬ 
odic Sobolev spaces and periodic Nikol’skij-Besov spaces with dominating mixed 
smoothness in [ 3 . They obtained c(s, a) = (s — l)/2 for order 2 digital nets in 
the former spaces for instance. Although the result is best possible, there are 
restrictions that only periodic function spaces are taken into account and that 
the smoothness parameter a, which equals r in their notation and is considered 
to be a positive real number, should be less than 2. Thus, the question arises 
whether higher order digital nets can achieve the best possible convergence of 
the integration error in non-periodic function spaces of a > 2. In this paper we 
give an affirmative answer to this question. 

To state the main result of this paper, we introduce some notation here. Let 
N be the set of positive integers and No := N U {0}. Let 6 be a prime and F;, 
the finite held with b elements, which is identihed with the set { 0 , 1 ,..., 6 — 1 } 
equipped with addition and multiplication modulo b. For x € [0,1), its 6 -adic 
expansion x = with is understood to be unique in the 

sense that inhnitely many of the ^i’s are different from 6 — 1. The operator 0 


2 


denotes digitwise addition modulo b, that is, for x = ^ = 

E“i € [ 0 > 1 ). we define 

OO 

x®x':= ^ ^ with rji = + (mod h). 

i=l 

Note that x (B x' is not always defined via its unique &-adic expansion and 
even may equal 1 ^ [0,1). Such an instance is given by setting b = 2, x = 
2~^ + 2~^ + 2“^ + • • • and x' = 2“^ + 2~^ + 2“® + • • •. However, if either x or 
x' can be written in a finite 6 -adic expansion, this situation never occurs. 

Moreover, let H be a normed function space with norm Il'Hy. The worst-case 
error of QMC integration over P in V is defined as 

e-”(H;P):= sup |/(/; P) - J(/)|. 

/ev 

ll/llv<l 

For (T G [0,1)®, we write P (B cr := {x (B cr : x G P}, where © is applied 
componentwise. Since we shall only consider a point set P whose each element 
X can be written in finite 5-adic expansions in this paper, a: © cr is always 
defined via unique 6 -adic expansions. Then the root mean square (RMS) worst- 
case error of QMC integration over P © cr in H with respect to a randomly 
chosen cr G [0,1)® is dehned as 

grms-wor(^.p) _ I f (e™''(H; P © cr))^ dcT 

Now the main result of this paper is given as follows. 

Theorem 1. For a, /?, m G N and t G Nq with a>2,j3>2a and 0 < t < ^m, 
let P be an order fd digital [t,m, s)-net owerFf,. Let 7 = (7u)uc{i,...,s} « set 

of the weights. The RMS worst-case error of QMC integration over P ® cr in 
TLa,~f with respect to a randomly chosen cr G [0,1)® is bounded by 

e“”(P„.^;P)<^ ^ 7y^a©.M.«aogiV)(l“l-')/^ ( 1 ) 

where Ca,p,b,t,u > 0 for all ^ u C { 1 ,..., s} and = jPj = 6 "*. 

Note that the explicit form of Ca,p,b,t,u can be found later in Q. This result 
directly implies the following. 

Corollary 1. For a, /3, m G N and t G No with a > 2, (3 > 2a and 0 < f < (dm, 
let P be an order jd digital [t,m, s)-net overWb- Let 7 = (7u)ttc{i,...,s} « set 

of the weights. There exists a cr G [0,1)® such that the worst-case error of QMC 
integration over P © cr m TLa.-y is bounded by 

® <^) < ^ E Ar)(l«l-i)/ 2 . 

0/uC{l,...,s} 

Although the t-value and thus the constants Ca,g,b,t,u may depend on m, it 
was shown by Dick Hi] that for large m we can explicitly construct an order 
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j3 digital it, m, s)-net with its t-value independent of m, see also Remark [T] 
Therefore, the rate of convergence which we obtain in this paper is of order 
iV““(log unless 7 {i,...,s} = 0. This compares favorably with what was 

obtained by Baldeaux and Dick in [TJ Theorem 24], where they considered the 
case where P is an order a digital net over F{,, i.e., the case where the order 
of digital nets and the smoothness parameter coincide, and obtained a similar 
bound on the RMS worst-case error but with the exponent of the logarithmic 
term equal to sa/2. Our result shows that the exponent of the logarithmic 
term is actually independent of a. Here we note that the convergence of order 

iV -“( logiV )('’-^)/2 

in a similar function space has been proven by using the 
Frolov cubature rule in conjunction with periodization strategy, see for instance 
| 20 ) . which is not a QMC integration rule though. 

Moreover, from the results of [5], we can see that the above results in Theo- 
rem[T]and Corollary [T] are best possible. Now let P = {a^o, ■ • ■ C [0,1)® 

be an N element point set and w = {wg, ■ • ■, rcAr-i} be an arbitrary real tuple. 
The worst-case error of cubature rule with points in P and weights w is defined 
as 




sup 

/ev 

ll/llv<l 


N-l 

Wnf{Xn) - I{f) 


n—0 


A lower bound on P, to) in the so-called half-period cosine space of 

smoothness a for the case of product weights, i.e., weights of the form = 
rijGu 7i) • ■ • 5 7s > 0; was proven in [6j Theorem 4], and furthermore, it 

was shown in [SI Theorem 1] that the half-period cosine space is continuously 
embedded in the Sobolev space 'Ha,-Y which we consider in this paper. Combin¬ 
ing these two results, we immediately have the following. 

Proposition 1 . Let a > 2 be a positive integer and 71,... ,7s > 0. For u C 
{!,..., s}, let ju = rijeuTj where the empty product equals 1. For any N 
element point set P C [0,1)'* and any real tuple w we have 






where Ca.-y.s is positive and independent of P and w. 

This implies that the exponent c(s, a) cannot be less than (s — l)/2 in 'Ha,'y 
for the case of product weights with 71 ,... , 7 s > 0. Since Corollary |T] shows 
the existence of point sets which achieve exactly this order, our result is best 
possible. However, as our result is again an existence result and thus is not fully 
constructive, it is interesting to study an explicit construction of deterministic 
point sets which achieve the best possible convergence of the worst-case error 
in Pa.-y- We leave it open for future work to address. 

In the next section, we shall introduce the necessary background and nota¬ 
tion such as weighted Sobolev spaces of smoothness a and higher order digital 
nets. In Section [31 we shall prove Theorem [H i.e., an upper bound on the RMS 
worst-case error of randomly digitally shifted order P digital nets in 
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2 Preliminaries 


2.1 Weighted Sobolev spaces 

First let us consider the one-dimensional unweighted case. The Sobolev space 
which we consider is given by 

Ha :={/: [ 0 , 1 )^K|/W: 

absolutely continuous for r = 0,..., a — 1, G T^[0, 1)|, 

where denotes the r-th derivative of /. As in pri Section 10.2] this 
space is indeed a reproducing kernel Hilbert space with the reproducing ker¬ 
nel ICa ■ [0,1) X [0,1) —;■ R and the inner product (•, ■)a given as follows: 

r Br{x)Br{y) , ,^^+i B 2 a{\x - y\) 

(^ 1)2 ( 2 a)! 

for x,y G [0,1), where Br denotes the Bernoulli polynomial of degree r, and 


{f,9)c,=‘^ ( f^"'\x)dx [ g^'^'>{x)dx+ [ f^°‘\x)g^^'>{x)dx, 

^^0 -^0 -^0 -^0 

for f,g € Ha. 

Let us move on to the s-dimensional weighted case. In the following we 
write {1 : n} := {l,...,n} for n G N. Let 7 = ( 7 u)uc{i:s} be a set of non¬ 
negative real numbers which are called weights. Note that the weights moderate 
the importance of different variables or groups of variables in function spaces 
and play an important role in the study of tractability m- However, such 
an investigation is out of the scope of this paper since we are interested in 
showing the optimal exponent of logiV term in the error bound. We consider 
the weighted function space for the sake of completeness. Moreover, we shall 
use the following notation: For u C {1 : s} and x G [0, 1)®, let Xy = {xj)j^y. For 
C u C {1 : s} and ry\y = irj)j^u\v, i'>^u\v,Oiv,0) denotes the s-dimensional 
vector whose j-th component is Vj if j G u \ u, a if j G v, and 0 otherwise. Now 
the weighted Sobolev space Ha ,7 which we consider is the reproducing kernel 
Hilbert space whose reproducing kernel /Ca,^: [0,1)'* x [0,1)'* —>■ R and inner 
product (•, •)a .7 are given as follows [T]: 


Aa,7 : y') 


hC{1:s} j£u 



Hr )-^r iVj ) 

(Hp 


\a+l B2a{\Xj yj\) 

^ ’ ( 2 a)! 


for X = (xi,..., Xs), y = (?/i,..., t/s) G [0,1)**, where the empty product always 
equals 1 , and 


(/, 5 )a, 7 = E E 

uC{l:s} tiCu r„\„e{l:a-l}l“\«' 


/ 

./[ 0 . 1 )l”l 




jir 


,\v, 


, 0 ) 


{x) dx 


{1:s}\d 
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X 


for f,gG 'Ha,-y, where for u C {1 : s} such that 7 ii = 0 we assume 


(/. 









X ( / g(’’“\’”“’'’°^(a:)da;{i.s}\^ ) da;„ = 0. 

\J[0,1)‘-M J 

Note that an integral and sum over the empty set is the identity operator and 
we formally set 0/0 := 0. 

2.2 Higher order digital nets 

Here we start with the general digital construction scheme of point sets as 
introduced by Niederreiter [Si- 


Definition 1. For m,n,s G N, let Ci ,... ,Cs G 0 < h < b"^ be an 

integer with its b-adic expansion h = ■ For 1 < j < s, let us consider 



where , ^n,h.j are given by 

J ^2,h,j J ■ ■ ■ 1 ^n,h,j) = Cj (^0 1 j ■ ■ ■ i —1) • 


The set P = {xq; ®i, ■ • ■, Xbm_i} c [0,1)® with Xh = {xh,i, ■ • ■, Xh,s) is called a 
digital net over F;, (with generating matrices Ci,..., Cg)- 

The dual net of P, denoted by is defined as follows. 

Definition 2. For to, n, s G N and Ci,... ,Cs G F^ let P be a digital net 
over Fh with generating matrices Ci,..., Cg. The dual net of P is defined as 


:= {fc = (fci,..., ks) €Wo:Cjki(B---®Cjks=0& F^} 


where we set k := (kq, ..., k„_i) for /c G Nq with its b-adic expansion k = 
kq + Kib + ■ ■ ■, which is actually a finite expansion. 

For a G N, we define a metric function as follows. 

Definition 3 . Let a G N. For fc G N with its b-adic expansion k = + 

K 2 b^^~^ + ■ ■ ■ + Kyb^'^~^ such that G {1,...,6—1} and ci > C 2 > • • • > 

Cy > 0. Then we define 



and gaiO) ■= 0. For k = (fci,..., kg) G Nq, we define 


S 
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Note that the above definition was originally given in [HIIS] for the case a = 1 
and in [Sill] for a > 2. We simply call /i^ the Dick metric function for any 
a > 1 throughout this paper. Now we define the minimum Dick metric of a 
digital net, which shall play a critical role in the subsequent analysis. 

Definition 4. Let P be a digital net over Ff, and its dual net. For a S N, 
the minimum Dick metric of P is defined as 

Sa(P) := min Ua(k). 

feGP-L\{0} 

Now we give the definition of higher order digital nets. 

Definition 5. For m.,n,a,s G N with n > am, let P be a digital net over Fj, 
with generating matrices Ci,... ,Cs G For 1 < i < n and 1 < J < s, we 

denote by Cij G FjJ* the i-th row vector of Cj. Let t be an integer with 0 < t < 
am which satisfies the following condition: For all 1 < ij^y. < ... < ij^i < n 
with 

s min(a,Vj) 

<am-t, 

3=1 1=1 

the vectors ,i, ■ ■ ■, ^ 1 , 1 , 1 , • • ■ 1 ..., are linearly independent over 

Fft. Then we call P an order a digital {t,m, s)-net over Ff,. 

The following property of order a digital (t, m, s)-nets directly follows from the 
linear independence of the rows of the generating matrices, that is, for any order 
a digital (t, m, s)-net P over Ff,, we have 

Sa{P) > ctm — t. 


Moreover, the following lemma is an obvious adaptation of the result shown in 
[SI Theorem 3.3] and [H Theorem 4.10], which states that any order a digital 
net is also an order a' digital net as long as 1 < a' < a. 

Lemma 1. For a € N, let P be an order a digital (t,m, s)-net over F;, with 
some integer 0 < t < am. Then, for any a' S N with I < a' < a, P is also an 
order a' digital {ta' ,m, s)-net overWb with t a' = \ta'/a]. 

Dick [SI S] proposed the following digit interlacing composition to obtain 
explicit construction of higher order digital nets over Fh: For m,s,a G N, let 
Q C [0,1)“® be a digital net over F;, with generating matrices Ci,..., Cas G 
For 1 < I < m and 1 < j < as, we denote by c^j- the i-th row vector 
of Cj. We now construct a digital net P C [0,1)® over Ff, with generating 
matrices Di,..., Dg G such that the (a(h — 1) + i)-th row vector of Dj 

equals Cb^a(j-i)+i for all 1 < h < m, 1 < i < a and 1 < j < s. Regarding this 
construction algorithm, we have the following, see for instance [2j Corollary 3.4]. 

Lemma 2. Let Q be an order 1 digital ft', m, as)-net over Ff, with 0 < t' < m. 
Then a digital net P constructed as above is an order a digital (t,m, s)-net over 
Ff, with 


t = a min 


m, t' + 




( 2 ) 
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Thus in order to obtain an order a digital (t, m, s)-net with small t-value, we 
need an order 1 digital (t',m,as)-net with small t'-value. Here we recall that 
there have been many explicit constructions of order 1 digital sequences (defined 
below) over F;, for arbitrary dimension proposed in the literature, so that we 
can construct order 1 digital as)-nets with small t'-value. 

Definition 6 . Let Ci,... ,Cs G 6e N x N matrices over Ff,. For Cj = 

(cj,fc,j)fc,/eN we assume that there exists a function /tT : N —?► N such that Cj^k,i = 
0 when k > K{1). Let h be a non-negative integer with its b-adic expansion 
h = rjiF for some a £ N. For 1 < j < s, let us consider 




where ■ ■ ■ are given by 

— Clj(r]Q,rii, ... ,r]a— 1 , 0 , 0 ,...) . 

The sequence S = {xo,Xi,...) with Xh = {xh.i,. ■. ,Xh,s) is called an digital 
sequence over F^ (with generating matrices Ci,..., Cg). 

Moreover, let t be a non-negative integer. For m £ N, let Cj^mxm be 
the the left upper m x m sub-matrix of Cj. If for all m > t the matrices 
Cl: ,771 X m 7 ■ • ■ 5 Cs.TnxTn generate an order 1 dijgital {t,m, s) net over F^, we call S 
an order 1 digital it, s)-sequence ouerFf,. 

Remark 1 . As mentioned above, there are many explicit constructions of order 
1 digital sequences over F;,, see for instance We refer to 0 

Chapter 6] for more information on this topic. For a, s £ N, let S he an order 1 
digital ft', as)-sequence over F^ with generating matrices Ci,..., Cas for some 
non-negative integer t'. Now let us define mo := t' + [s(a — 1)/2J. When 
m > mo, by using the result of Lemma\^ we see that the digital net P C [0,1)® 
constructed by the digit interlacing composition based on Ci^rnxm, ■ ■ ■, Cas.mxm 
becomes an order a digital (amo,m, s)-net. Here the value amo does not depend 
on m. 


3 Proof of Theorem [T] 

Throughout this section, let P be an order (3 digital net over Ff, for /3 £ N. Here 
we prove Theorem[Tl i.e., an upper bound on the RMS worst-case error of QMC 
integration over P 0 cr in Ha.-y with respect to a randomly chosen a £ [ 0 , 1 )® 
when (3 > 2a. 

3.1 Interpolation of Dick metric functions 

In this subsection, we discuss an interpolation property of Dick metric functions, 
which shall become a crucial tool in the proof of an upper bound on the RMS 
worst-case error. 

Lemma 3. Let a,/? £ N with 1 < a < fi. For any k £ Ng, it holds that 

/.N (3 — a ,,, 

Ma(fc) > 
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For convenience, we shall in what follows write 


^ck/3 


a — 1 


and i?o ,/3 


P — a 

J^' 


Proof. Since fj,a{k) = ^J■a{kj) for any a € N and k = (fci, ... ,ks) € Nq, it 

suffices to prove that the inequality 


^ ^aj3 “t" -Sck/3(^); 

holds for any k gNq. As the result is trivial for k = 0, we only consider the case 
A: > 1 in the following. Let us denote the 5-adic expansion oi khy k = Kib'^^~^ + 

H-h for some u > 1 such that ki, ..., G {1 ,..., 6 — 1} and 

Cl > • • • > c„ > 0. When j5 > v, we write c„+i = c „+2 = • • • = c^ = 0. Then, 
we have 

a. OL 

y.a{k) = ^Ci > Cl + ^Cq, = /ri(fc) + (a - l)Ca, 

i=l i=2 

as well as 


0 a P 

Hisik) ='^Ci ^ Ca = fJ.a{k) + if) - a)Ca. 

2=1 2=1 2 =Q ;+1 

By using the above two inequalities, we obtain 

Hj3{k) - Haik) ^ ^ ^ ^j,a{k) - ni{k) 

o S Ca S , 

p — a a — 1 

from which we can easily see that the result follows. □ 


Remark 2. In the proof of the upper bound on the RMS worst-case error, which 
shall be given in the next subsection, the condition Bap > 1/2 is necessary. This 
condition can be satisfied if and only if ft > 2a. This is why we assume ft > 2a 
in Theorem]^ and Corollary]^ 


3.2 Upper bound on the RMS worst-case error 

Finally, we prove an upper bound on the RMS worst-case error of QMC inte¬ 
gration over P 0 cr in TLa.-y with respect to a randomly chosen cr G [0,1)'*. The 
following lemma stems from the proof of [TJ Theorem 30]. 

Lemma 4. Let P be a digital net over and P^ its dual net. For u C {1 : s}, 
we write P/p = G : (fc„,0) G P"*"}- The mean square worst-case error 
of QMC integration over P (B cr in TLa.-y with respect to a randomly chosen 
cr G [0,1)® is bounded by 

07 ^: 2 iC{l:s} kuGP.jjf 


where we simply write Paiku) := Ta(,kj) for 0 7 ^ it C {1 : s} and k^ G 

and Da. b > 0 depends only on a and b and is given by 


Da,b = max 

’ l< 2 ;<a 



{Cr,b? 

b^C-v) 


‘2C2a,b I 

ip2(a-v) I ’ 
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with 


Ci.b = 


1 


and Crb = 


(l + l/6+l/(6(6+l)))"-2 


for T > 2. 


2 sin(r/&) ’ (2 sin(r/6))'^ 

In the subsequent analysis, we shall use the following inequality, see [TJ 
Lemma 13.24] for its proof. 

Lemma 5. For any real number b > 1 and any k,tQ G N, we have 

-k 


t + k — 1 
k-1 


< 5 -to ( ^0 + fc - 1 
k — 1 


t=to 

Now we are ready to prove Theorem [T] 


1 - 


1 


Proof of Theorem [IJ Using Lemmas [3] and SJ we have 




< 'Y^ Y^ ^-‘2'AccpSp{P)-2Bccpni{ku) 

0^LiC{l:s} 

07^«C{l:s} 


where 5p{P) is defined as in Definition 0] and we write 

WP^'^%P):= Y 

for 0 7 ^ u C {1 : s}. Since /3 > 2a, we have Ba^ > 1/2 as stated in Remark [2] 
In the following we focus on the term (P). Since pi(fcu) is an integer 

no less than both juj and Si{P) for any fc„ G we have 


OO 

wP^-^{P)= Y E 

?i=max{(5i(P).|«|} k,^ePf 

OO 

b-^Bo.ph ^ 

?i=max{(5i(P),|u|} 

Ul(k-u.)=h 
OO 

E E E 

/!,=max{(5i(P).|u|} k,^ePf 

\iu\i—k yjeu 



where we denote \lu\i = P,j^u h- innermost sum in the last expression, 

it is known from m Lemma IS.sJ^ that we have 

fO if |Z„|i <5i(P), 

^ l<i(6-l)H if<5i(P)<|l.|i<5i(P) + |w|, 

kuGPf [(5_ i)Dl5lD|i-(5i(P)+DI-i) if |Z„|i > 5i(P) + |m|. 

^Although [3 Lemma 13.8] only consider the case where P is a digital net over with 
generating matrices of size m x m, the proof still goes through even when P is a digital net 
over Ff, with generating matrices of size n x m as long as n > m. 
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Thus {P) can be bounded by 


5i(P) + |«|-l 

/i=max{5i(P),|ii|} i„eNl“l 

\K\l=h 
oo 

;i=5i(p)+|u| 

\l^\l = h 

5i(P) + |«|-l 


= ( 6 - 1)1 


u| 


A—max{(5i (P), |ii|} 


h — 1 
m| — 1 


_|_ 5-('5i(f’) + l“|-l) Y^ 

h=Si{P) + \u\ 


h-1 

Iu|-1 


(3) 


For the second sum in Q, we have 


h-1 

| m |-1 


h=Si{P) + \u\ 

= Y b-C^B^ff-^){h+\u\)fh+\u\-l 

h=Si{P) 


< 


< 


V H-1 

^-(2B„p-l)(5i(P) + |u|) (l- f,-{2B^f>-l)'^ 


-I“l 


1 


l“l 


((5i(F) + l)l“l-i 

pB^ls-l _ij 5(2B„3-1)5i(P) ^ 


where we used Lemma 0 in the first inequality as we have Ba /3 > 1/2 by the 
assumption /3 > 2a, and the second inequality stems from the inequality 

a.(P)^l l-n 

V |w|-i / 

For the first sum in JS]), we have 

5i(P) + |a|-l 

y- u-2B^phl n i 

^ llul - 1 

/z=max{(5i (P), liil} 


< "Y^ Ij-'iBafih 

/i—max{5i (P), |ii|} 


h-1 

Iu |-1 


< 


Y l^-2B^f,ih+\u\) fh + \u\ - 1\ 

/i—max{5i (P) —l^zl ,0} ^ ^ 

max{5i(p).|„|} f max{6i(P), \u\} - 1\ 

V |w|-i 


-\u\ 


where we used Lemma [5] again in the last inequality. 
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Now let us consider the case Si{P) > |m|. The first sum in ^ is bounded by 


^i(-P)+l“|-i 

/i—max{(5i (P), |ii|} 


h- 1 
u\ — 1, 


< 


< 


1 b2Ba.fi _ 1 


/ y^Ba^fi Y“l (<5^(P) _ l)l«|-l 


IpBafi&li.P) 


b^Bcfi _ 1 


b-^Ba^fi \ l“l (J^(p) + l)l«|- 


b'2Bo,fiSi{P) 


For the case Si{P) < |u|, the first sum in (I2|) is bounded by 

(5i(P) + |«|-l / ^ 1 \ 

b-^Ba.fihlB- M < 

JR. ( D^ U,l\ Vl^l 


/i=max{(5i (P) ,|n|} 


< 


lp'Ba0 — ^ 


b^Ba^fi \ l“l (5i(P) + l)l“|-l 


b2Ba,fiSi{P) 


Thus, regardless of whether Si{P) > |u| or 5i{P) < |w,|, we have the bound on 
the first sum in Q as 


5i(P) + |«|-l 
/i—max{(5i (P) ,|u|} 


h — 1 
mI — 1 


< 


b‘2-Bc.fi — 


b^B.,fi \ l“l 1)N|-1 


b2Ba,fiSi{P) 


Combining this result with the bound on the second sum, we have 


WyBa.fi ^p-^ < ( 6 - 1 ) 


|tt| 


yb^Bafi _ 1 

+ ;,-(<5i(P)+NI-i) 




b2Ba,fiSi{P) 


1 


l“l 


(5i(P) + 1)H 


-1 


lpBct0 — l — 1 y l){‘^Ba0 — l)6i{P) 


= G. 


((5i(P) + l)l“l-i 

a,P,b,u ^2BafiSj{P) ’ 


where we set 

Gc,p,b.u = {b-l)\^\ 

So far we have obtained 


yZBafi 
62Bo/3 _ 1 


l“l 


+ b 


1 


b'^Bafi _ b 


Nl 




0/ixC{l:s} 


< 


1 


b2Aa.fiSfi{P) + 2BafiS\{P) 


Y. 7«<IG„,;5P„((5i(P) + 1 )I“|-1. 


0 y^iiC{l:s} 
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Finally let us recall that P is an order /3 digital (t,m, s)-net over F;,. From 
this fact and Lemma [TJ we have 

Si{P)>m — ti and 5p{P)>l3m — t, 

where ti = Thus, we have 

2Aaf3Si3{P) + 2Bap5i{P) > 2Aap{l3m -t) + 2Bap{m - h) 

= {2f3Aap + 2Bap)m — 2Aapt — 2Bal3ti. 

In the above, it holds that 


2pAap + 2Bcti3 — 2 a. 

Since ti = 0 is best possible, we have 6i{P) < m+1. Therefore, we get 








Jj2am 

0/uC{l:s} 

L2Aa^t-\-2Bct^ti _ 

05^mC{1:s} 


62° 


which completes the proof by using the inequality ^ 

tti > 0 and by choosing 


, ,, , / o \(l«|-l)/2 

'^a,/3,b,u,t — 0 ^a,b ^a,l3,b,u I Jog 6 / 

for all 0 7 ^ u C {1 : s} such that the bound ((T|) holds. 


(4) 

□ 
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